Non-thermal origin of nonlinear transport across magnetically induced 
superconductor- metal-insulator transition 
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We have studied the effect of perpendicular magnetic fields and temperatures on the nonlinear 
electronic transport in amorphous Ta superconducting thin films. The films exhibit a magnetic field 
induced metallic behavior intervening the superconductor-insulator transition in the zero tempera- 
ture limit. We show that the nonlinear transport in the superconducting and metallic phase is of 
non-thermal origin and accompanies an extraordinarily long voltage response time. 
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In recent years, the suppression of superconductivity 
in two-dimensions (2D) by means of increasing disor- 
der (usually controlled by film thickness) or applying 
magnetic fields has been a focus of attention. Conven- 
tional treatments Q, |^ Q, |^ of electronic transport 
predict that in 2D the suppression of the superconduc- 
tivity leads to a direct superconductor-insulator tran- 
sition (SIT) in the limit of zero temperature (T = 0). 
This traditional view has been challenged by the obser- 
vation of magnetic field (B) induced metallic behavior in 
amorphous MoGe 6^ 1_, _8J and Ta thin films P|. The 
unexpected metallic behavior, intervening the i?-driven 
SIT, is characterized by a drop in resistance {p) fol- 
lowed by a saturation to a finite value as T — > 0. The 
metallic resistance can be orders of magnitude smaller 
than the normal state resistance (p„) implying that the 
metallic state exists as a separate phase rather than 
a point in the phase diagram. Despite many theoret- 
ical treatments |i[lS[lllIIlIIllllIIlIllII3, a 
consensus on the mechanism behind the metallic be- 
havior is yet to be reached. Proposed origins of the 
metallic behavior include bosonic interactions in the non- 
superconducting Dhase|ilollllj. contribution of fermionic 
quasiparticles to the conduction^, 0|, and quantum 
phase fluctuations 0, 01 . 

In a recent paper P| on the magnetically induced 
metallic behavior in Ta films, we have reported the non- 
linear voltage-current {I-V) characteristics that can be 
used to identify each phase. The superconducting phase 
is unique in having both a hysteretic I-V and an "im- 
measurably" small voltage response to currents below 
an apparent critical current Iq. The metallic phase can 
be identified by a differential resistance {dV/ dl) that in- 
creases with increasing /, whereas the insulating phase 
is identified by a dV/dl that decreases with increasing 
/. The contrasting nonlinear I-V in the metallic and 
insulating phase are shown in Fig. 1(a). 

The main purpose of this Letter is to report that the 
origin of the nonlinear transport, particularly in the su- 
perconducting and metallic phase, is not a simple reflec- 
tion of T-dependence of p via the unavoidable Joule heat- 
ing. We describe the effect of B and T on the nonlinear 



TABLE I: List of sample parameters: nominal film thickness, 
mean field Tc at i? = 0, normal state resistivity at 4.2 K, crit- 
ical magnetic field at which the resistance reaches 90% of the 
high field saturation value, and correlation length calculated 
from ^ = ^ where "l>o is the flux quantum. 



sample 


batch 


t{nm) 


Tc{K) 


Pn(!^/n) 


Be 


C 


Ta 1 


1 


3.5 


0.584 


1769 


0.72 


21 


Ta 2 


1 


5.0 


0.675 


1180 


0.88 


19 


Ta 3 


2 


5.7 


0.770 


1056 


0.9 


19 


Ta 4 


3 


5.0 


0.598 


770 






Ta 5 


4 


36 


0.995 


69 


2.0 
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transport on which this conclusion is based. We also 
present our studies on dynamic voltage response which 
reveal strikingly long voltage response times that accom- 
pany the nonlinear transport. 

Our samples are dc sputter deposited Ta thin films 
on Si substrates. The sputter chamber is baked at 
~ 110° C for several days, reaching a base pressure of 
10-8 Torr. The chamber and Ta source were cleaned by 
pre-sputtering for ^ 30 min at a rate of ~ Inm/s. Films 
are grown at a rate of ^ 0.05 nm/s at an Ar pressure of 
'-^ 4 mTorr. Using a rotatable substrate holder up to 12 
films, each with a different thickness, can be grown with- 
out breaking the vacuum. In order to facilitate four point 
measurements, the samples are patterned into a bridge 
(1 mm wide and 5 mm long) using a shadow mask. Even 
though there were noticeable batch to batch variations, 
the degree of disorder (evidenced by the values of p„) for 
films of the same batch increases monotonically with de- 
creasing film thickness. The superconducting properties 
of the films are characteristic of homogeneously disor- 
dered thin films and consistent with the results of 
x-ray structural investigations . The data presented in 
this paper are from 5 films grown in 4 batches. Parame- 
ters of the films are summarized in Table HJ 

The evolution of the I-V curves across the 
superconductor-metal boundary at 40 mK is shown in 
Fig. ^b) for sample Ta 1. The hysteretic I-V, unique 
to the superconducting phase, is indicated by the dashed 
lines. As / is increased, the superconductivity is abruptly 
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FIG. 1; (a) dV/dl vs. / at 20 mK across the metal-insulator 
boundary for a Ta film with Tc = 0.23 K, adapted from Ref.0|. 
(b) Current biased I-V curves of sample Ta 1 at 40 mK and 
the indicated B. (c) Our method of I-V measurements is illus- 
trated, (d) The T-dependence of p for sample Ta 2 at _B = - 
1.0 T with 0.1 T interval and 3.0 T, measured at 7 Hz with a 
current amplitude of 1 uA. The solid lines are to indicate the 
superconducting phase in the low T limit, the dashed lines the 
metallic phase, and the dotted lines insulating phase, (e) I-V 
curves in log-log scale for sample Ta 3. Filled triangles (open 
circles) are for current increasing (decreasing) branch. The 
data density is reduced to make individual symbols visible. 
The arrow marks the critical current Ic and voltage Vc. 



quenched at a well-defined critical current Ic- As I is de- 
creased from above Ic, the superconductivity suddenly 
appears at a different current I'c < Ic- The hysteresis be- 
comes smaller with increasing B, and vanishes near 0.1 
T as the system is driven into the metallic phase [solid 
lines in Fig. db)]. Typical T-dependence of p at various 
B is shown in Fig. ^d) for another sample Ta 2. In this 
sample, which is less disordered than Ta 1, the supercon- 
ducting phase extends up to ^ 0.2 T [solid hues in Fig. 
n^d)] and the metallic behavior is observed at higher B 
up to ~ 0.9 T (dashed hues). Hereafter, "superconduct- 
ing regime" refers to the transport regime where I-V is 
hysteretic and "metallic regime" to the regime with non- 
linear (and reversible) I-V with increasing dV/dl with 
increasing /. 

All our I-V curves are constructed by plotting the 
steady state voltage at each bias current that is changed 
in small discreet steps as illustrated in Fig. dc). In or- 
der to ensure that the steady state is reached at each 
step, the voltage is monitored every 50 ms for up to 55 
s while the current is kept constant. The magnitude of 
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FIG. 2: (a) I-V curves of sample Ta 4 at B = and at the 
indicated T. Open (filled) circles are for the current increasing 
(decreasing) branch, (b) T-dependence of p for sample Ta 4 
at B = 0. The arrow is to indicate the temperature T* where 
the hysteresis vanishes, (c) The I-V of the sample Ta 5 at B 
= 0. The temperature of each trace is, from the top, 1.000, 
0.994, 0.992, 0.990, 0.988, 0.984, and 0.976 K. (d) The T- 
dependence of p for sample Ta 5 at B = 0. 



the voltage jump at Ic (or /^), which could be as large as 
several orders of magnitude, was almost independent of 
the current step size in the range 5 - 100 nA. Even with 
our smallest steps of 5 nA, no steady state with a voltage 
within the range covered by the jump was observed. 

Our investigations on how B and T influence the non- 
linear transport indicate that these quantities play simi- 
lar roles. The main effect of increasing T is to lower the 
superconductor-metal "critical" field i?*™; the i3-driven 
evolution of the I-V curves at an elevated T [Fig. ^e)] 
remains qualitatively the same as that in the low T limit 
[Fig. nib)]. More importantly, the evolution of the I-V 
curves as a function of T [shown in Fig. \^a) and (c)] 
is strikingly similar to that caused by B [shown in Fig. 
mb) and (e)]. 

The field B^™ decreases with increasing T and reaches 
zero at a well-defined temperature T*, which is close to 
Tc as shown by the arrows in Fig. |2Ib) and (d). This, 
together with the observations described above, means 
that i?^™ is a well-defined line in B-T plane separat- 
ing the superconducting and metallic regime. We point 
out that, in terms of nonlinear transport, the electronic 
properties at B > S^™ in the low T limit, where the un- 
expected metallic behavior intervening SIT is observed, 
are indistinguishable from those at high temperatures, 
for example at T > T* and S « 0. 

An important finding from our investigations is that 
the voltage jump in the superconducting regime is of 
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FIG. 3: (a) The B-dependence of the critical power, Pc=^cV'c, 
for sample Ta 3 at 0.5 K. 7c and Vc are marked by an arrow 
in Fig. 1(e). The dashed line is to guide an eye. (b) vs. B 
plot for sample Ta 5 at 0.970 K. 



non-thermal origin. If Joule heating is significant, the 
electron temperature would be determined by the bal- 
ance of the Joule heating power and the heat drain rate 
to the stage where the sample is thermally anchored. A 
large heating power at a high bias current could make 
Te substantially higher than the stage temperature. If 
Te reaches near where p sharply rises with T, an in- 
crease in Te could cause an increase in the heating power, 
which in turn causes a further increase in Tg. Such a 
positive self-feedback would make Te unstable and run 
away beyond Tc, resulting in a sudden quenching of the 
superconductivity appearing as a voltage jump. This sce- 
nario can be tested by applying weak B. The magnetic 
fields lower while the net thermal conductance be- 
tween the sample and the stage would remain almost un- 
affected. Therefore, in the heating scenario the critical 
power, Pc — IcVc where /c and Vc are the current and 
voltage at the onset of the voltage jump on the current 
increasing branch [marked by an arrow in Fig. ^e)], is 
expected to be weakly decreasing function of B. How- 
ever, as shown in Fig. Ola) and (b), is found to in- 
crease by an order of magnitude or more under weak B. 
This clearly demonstrates that the voltage jump in the 
superconducting regime has a non-thermal origin. We 
note that Vc is the highest steady state voltage in the su- 
perconducting state before the current-induced sudden 
quenching of the superconductivity. In repeated runs af- 
ter thermal cyclings to above 10 K, the value of Vc was 
reproducible within several percent even with different 
current step sizes in the range 5 - 100 nA. 

The 1-V curves in Fig. I1I2I clearlv show that the dis- 
continuity in 1-V in the superconducting regime evolves 
into the point of the largest slope in the continuous 1-V 
in the metallic regime. This strongly suggests that the 
sudden quenching of the superconductivity at Ic and the 
nonlinear transport in the metallic regime are caused by 
the same mechanism, which has been shown above to be 
of non-thermal origin. 

Now we turn to the discussion on the dynamic voltage 
response. The dynamic voltage response was studied by 
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FIG. 4; (a) An example of V-t trace at / = 24.5 A. The solid 
line is a fit to the exponential function shown. The dashed 
line is a fit to a logarithmic i-dependence for the first 3 s. (b) 
t vs. / plot at the indicated B and T. Filled (open) circles are 
for current increasing (decreasing) branch. The dashed lines 
are to guide an eye. Inset: the slope of 1-V curves in log-log 
scale at the same B and T. (c) r-J plot at _B = 0. The tem- 
perature of each trace is, from the top, 0.984, 0.988, 0.990, 
0.992, 0.994, and 1.000 K. Each trace is vertically shifted suc- 
cessively. For the bottom four traces, only current increasing 
branch is shown. Inset: d(log y)/d(log /) is plotted at the 
same T and B. Each trace is vertically shifted successively. 



analyzing the voltage-time (T^-t) traces to determine how 
fast the steady state is reached at each bias current. As 
shown by the solid line in Fig. ^a), the i-dependence 
is well described by an exponential function. The pa- 
rameters Vo and V\ are determined from the measured 
steady state voltage and V(t = 0) which is the steady 
state voltage at the previous bias current. The param- 
eter T is defined as the voltage response time constant, 
and obtained from a least squared error fitting procedure. 

Figure Elb) shows the t-1 plots at three different B 
across the superconductor-metal boundary. In the su- 
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perconducting regime (B=100 G), the time constant ex- 
hibits a hysteretic and diverging behavior with approach- 
ing /c from below (7^ from above) where the I-V is dis- 
continuous. In the metaUic regime (150 and 200 G) where 
I-V is continuous and reversible, the r-traces are also 
continuous and reversible. However, a prominent peak 
structure is evident in the data. Interestingly, the peak 
or diverging feature in r almost exactly coincides with 
the nonlinear transport. The inset shows the slopes of 
I-V curves in log-log scale. At current where r is large 
the transport is nonlinear [d{\ogV)/d{\ogI) > 1], and 
at current where r is almost zero the transport is linear 
[c?(log V) /c?(log /) = !]. Qualitatively the same r behav- 
ior is observed when the superconductor-metal boundary 
is crossed by increasing T. This is shown in Fig. ^c). 
Note that the traces in Fig. 2fc) and the inset are suc- 
cessively shifted vertically. 

It is surprising to find that the time constant in the 
peak or diverging region is as long as several seconds. 
We emphasize that the long time constant is not due to 
instrumentation. This is best demonstrated by the sys- 
tematics of the data. Changing B or T systematically 
shifts the peak or divergence while outside the narrow 
peak or diverging region the time constants remain al- 
most zero. The time constants were measured to be the 
same within the scatter of the data for current steps of 5 
nA (not shown) and 100 nA (shown). 

In the past, hysteresis accompanying a long response 
time has been studied in the context of irreversible 
magnetic properties in type II superconductors |l9j . 
Large magnetic relaxation rates observed in magnetiza- 
tion measurements [2(Tj | are believed to arise from thermal 
activation of magnetic flux lines out of pinning sites. The 
depinning process leads to a redistribution of flux lines 
causing a change in magnetization with time. Although 
our hysteresis with a long response time is observed in 
transport, not magnetization measurements, it still may 
be possible to understand in terms of pinning-depinning 
of vortices. A pinning-depinning transition arises from 
the competition between disorder-induced pinning force 
and Lorentz driving force due to the bias current. Un- 
der such a competition, the vortex motion is analogous to 
the flow of sand grains in a sand pile , where the com- 
petition is between the jamming due to the granularity 
of the system and gravitational force. Indeed, hystere- 
sis (22] and slow relaxation rates ^2^ have been observed 
in granular flow under mechanical vibrations. Logarith- 
mic time dependence has been observed in measurements 
of relaxation of magnetization llij and experiments on 
granular flow [2^. The dashed line in Fig. Efa) is the fit 
of our data to a logarithmic time dependence, and de- 
scribes the data well for t < 3 s. Nevertheless, over the 
entire range of the data, the exponential function (solid 
line) fits our data better. 

Finally, we briefly discuss an interesting implication 
of the unusually long response time accompanying the 



nonlinear transport at i? = and x > T* . Nonlin- 
ear transport of 2D superconductors at B = is usu- 
ally understood in the framework of Kosterlitz-Thouless 
(KT) theory where the superconducting transition 
corresponds to a thermodynamic instability of vortex- 
antivortex pairs in 2D. In this picture, current-induced 
dissociation of vortex pairs in the superconducting phase 
is expected to lead to nonlinear transport in the fashion, 
V (X I^ with a > 3 |2^. It has been argued |2^ that in a 
real system finite size effect can induce free vortices alter- 
ing the power law I-V. The resulting I-V curves obtained 
in numerical simulations [2^ show a pronounced peak 
structure in c?(logI^)/d(log/) resembling those shown in 
the inset of Fig. ^c). However, the voltage response time 
which is observed to be as long as several seconds, is too 
long to be reasonable with the KT framework where the 
vortex dissociation is expected to occur in the time scale 
of quasiparticle scattering, typically ~ lO^^s |23|. At 
present, whether dynamics of KT vortices in the pres- 
ence of disorder can have such long time scales is not 
clear, and the understanding of the nonlinear transport 
requires the development of a theoretical framework for 
nonequilibrium dynamics. 

In summary, we have shown the magnetic field and 
temperature driven evolution of phase-identifying non- 
linear I-V characteristics of Ta thin films that exhibit 
an unexpected metallic phase intervening SIT in the low 
T limit. Our observations indicate that a non-thermal 
mechanism is behind both the nonlinear transport in the 
metallic phase and sudden current-induced quenching of 
the superconductivity in the superconducting phase. Our 
dynamic voltage response studies suggest a possible link 
of the metallic behavior to the dynamics of vortices in 
the presence of disorder. 
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